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At a telephone exchange, calls appear before a single trunk line in ac- 
cordance with a Poisson process of density X. // the trunk line is busy, calls 
are delayed. The call holding times are identically distributed, mutually 
independent, positive random variables with distribution function H(x). In 
this paper the distribution function of the delay and its moments are given 
for a stationary process and for three orders of service: (i) order of arrival, 
(ii) random order, and (Hi) reverse order of arrival. 

I. INTRODUCTION 

Let us suppose that in the time interval (0, oo ) calls appear before a 
single trunk line at times t\ , r% , ■ ■ • , r„ , • • • where the interarrival 
times t„ — t„_i (n = 1, 2, • • • ; r = 0) are identically distributed, 
mutually independent random variables with the distribution function 

r . . (l - <T Xx iix ^ 0, , 1X 

F(X) =\0 ifs<0, (1) 

that is, the input is a Poisson process of density X. If an incoming call 
finds the line free, a connection is realized instantaneously. If the line 
is busy, the call is delayed and waits for service as long as necessary (no 
defection). The holding times are identically distributed, mutually in- 
dependent, positive random variables with distribution function II (x) 
and independent of the input process. Such a service system can be 
characterized by the symbol [F(x),II(x),l] provided that the. order of 
sendee is specified. In this paper three orders of service are considered: 
(i) order of arrival (first come-first served), (ii) random order (every 
waiting call, independently of the others, and of its past delay, has the 
same probability of being chosen for service), and (Hi) reverse order of 
arrival (last come-first served). 

487 



488 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1963 

We are interested in finding the distribution function of the delay for 
a stationary process and for the three orders of service. We shall prove 
that if Xa < 1, where a is the average holding time, then there is a unique 
stationary process. 

Throughout this paper we shall use the notation 

*(. s ) = r e - sx dH(x) (9?(s) ^ 0) (2) 

and 

a k m r x k dll(x) (k = 0,1,2, ••■). (3) 

•'o 

In particular, a { = a is the average holding time. 

II. THE STATIONARY PROCESS 

Let us denote by £„ the queue size at time / = r n — 0, i.e., the nth 
incoming call finds £„ calls (either waiting or being served) in the system. 
Denote by x« the time needed to complete the current service (if any) 
at time t - r„ - 0. If £„ = then x» = 0. The vector sequence 

(£„ ,X»), n = 1, 2, ••• , 

is a Markovian stochastic sequence and has the same stochastic behavior 
for each order of service. We shall prove that if Xa < 1 then there exists 
a unique stationary distribution, whereas if Xa ^ 1 then a stationary dis- 
tribution does not exist. If Xa < 1 and ($i , xO has the stationary dis- 
tribution, then every (£„ , x-0 has the same distribution as the initial 
distribution. For the stationary process, let us introduce the following 
notation 

Plfc. =j] = Pi Ci-0,1, •••) (4) 

P!xn ^ x, £„ = j) = Pj(x) (x ^ 0, j = 1, 2, • • • ) (5) 
and 

n y (s) = f e-"dPj(x) (9t(«) fc <W - 1,2, ■••). (6) 

We shall prove the following theorem, due to D. M. G. Wishart: 1 
Theorem 1: If Xa < 1, then the stochastic sequence (£„ , Xn), « = 1, 2, 

• • • , has a unique stationary distribution which is given by 1\ = 1 — Xa 

and 
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u(s, z ) = Eii,( s y 

(7) 



= (1 - \a)\z(l - z) ( *(s) - *(X(1 - z)) 
z - *(X(1 - z)) \ s- X(l - z) 



foryi(s) ^ Omul \z\ ^ 1. 

Proof: If we express the distribution of (£ n+1 , x „+i) with the aid of 
the distribution of (fc, , x»), and assume that both (f«+i , %n+i) and 
(£n , Xn) have the same stationary distribution, and if we form Laplace- 
Sticltjes transforms, then we obtain that P and II _,•($) (j = 1,2, • • •) 
must satisfy the following system of linear equations: 

P = Po*(x) + En ( .(A)[ff\)f, 



(8) 



S — A (_ t=l 

N( s ) = MlIj-i(X) ~ IWs)] 
s — X 

+ xt>( y - * M] i «.) W xr 

for j = 2, 3, • • • and 

Po+ ZlIXO) = 1. (9) 

To prove (8) we use the following two facts: First, the probability that 
during a holding time no call arrives is given by 



*(X) = f e^dHix). 



Second, let p and 6 be mutually independent, positive random variables 
with distribution functions Pjp fS x) = P(x), and P{0 ^ x] = F(x) 
denned by (1). Write A = \d < P ). Then 



P{A\E\c-" (p - e) 
where 



..,, _ f" f e -<-*-» dydPM = Mn(x) - n(.)i 

n(s) = f e~ sx dP(x). 
•'o 
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Forming generating functions in (8) we obtain 
[s - X(l - z)]U(s,z) 

= wkm) + m^x) -*(»)! {ft + p( ^j:f (x 4 (10) 

If s = \(1 - z) in (10) then we get 

t/(M) + WX) -*(x(i -*))] {ft + y(x '*g; I f-'V o. (ID 

The comparison of (10) and (11) gives 

„. Xg[^(g ) -tt(X(l - z))\ Tja x n2) 

U(s,z) = [s _ x(1 _ z)]mx) _ nx{1 _ m UM. (12) 

By the first equation of (8), J7(X,¥(X)) = Po[l - *(*)], and if we put 
this into (11) we get 



Po(l - z)[*(X) -¥(X(1 -g))1 
Thus by (12) and (13) 

U(S ' Z) = \z - ¥(X(1 -V))][8 - X(l - «)J ■ 



„,- v Pod - z)mv -*ixu-zj;j na 

^'^ = «-*(x(i-«)) • (13 



xPozd - «)[*(«) - ^(x(i -- g))] (14) 



Since by (9) P + f/(0,l) = 1, it follows from (14) that P = 1 - Xa. 
Thus if Xa ^ 1, then the assumption that a stationary distribution 
exists leads to a contradiction, i.e., a stationary distribution cannot exist 
if Xa ^ 1. If Xa < 1, then there exists one and only one stationary 
distribution which is given by P„ = 1 - Xa and by (14). This proves 

(7). 
Remark. From (7) we obtain by inversion that for x ^ 

ZPyW 

(15) 
= (1 - Xa)Xz(l - z) r e -xa-.>« [H ( M + x) _ H ( u )] dUm 
z -¥(X(1 - z)) h 

Hence for x ^ 

EP,(x) = X f [1 - H(u)]du. (16) 

Accordingly if (£„ , Xn), w = 1, 2, • • •, is a stationary sequence, then we 
have for x ^ that 
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P{xn £*|&fclj=i/"[l- H{u)\ du, (17) 

i.e., if an incoming call finds the line busy, then the distribution function 
of the time needed to complete the current service is given by 

- / [1 - H(u)] du if x ^ 0, 

a J ' (18 ) 



H*(x) = 



if x < 0. 



Finally we also remark that the stationary distribution of £„ , n = 1,2, 
• • • , is given by the following generating function 

This follows from (7), because U(z) = P + C/"(0,2). 

III. THE DISTRIBUTION FUNCTION OF THE LENGTH OF A BUSY PERIOD 

A busy period is defined as a time interval during which the line is 
continuously busy. The stochastic law of a busy period is obviously 
independent of the order of sendee. Every busy period (except the 
initial one, if the line is busy at time t = 0) independently of the others 
has the same stochastic law. Denote by G(x) the probability that the 
length of a busy period (other than the initial one, if any) is ^x and 
define 

7 (s) = f e~ 8X dG(x) («R( S ) ^ 0). (20) 

In Ref. 2 it is proved that y(s) is the root with smallest absolute value 
in 2 of the equation 

z - *(s + X(l - z)). (21) 

By Lagrange's expansion (cf. Ref. 3, p. 132) we obtain that 

7(s) = E ^ f e-^'x-UHM (22) 

where H n (x) denotes the nth iterated convolution of H(x) with itself. 
From (22) it follows by inversion that 

G(x) = E^t- 1 I' e^u-'dHniu). (23) 

n=l Til ''O 

If \a ^ 1 then G( oo ) = 1, whereas if \a > 1 then G( » ) < 1. 
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In the case of Xa ^ 1 the rth moment of G(x) is defined by 
T r = r x r dG(x) (r = 0,1,2, •••). 



(24) 



If Xa < 1 and tt T is finite then T , r, , ■ • • , IV are also finite and we 
have r = 1, Ti = a/(l - Xa), and 



n+1 £tn!(l - Xa)»+' +1 
for n = I, 2, • • •, where 

v n! oi2 h az h ■ ■ • 

Yn * ~ ,•,+/,+"■+/.=> *l*l • • • i.!(2l)'»(3!)'« .-.((»+ 1) I)*. 



r„., 



OTn+l 



(25) 



rrr • (20) 



If, in particular, H(x) = 1 - e xla (x ^ 0), then a r = r\a and 



Yn., = « 



E 



n\ 






! /n — l\ „+, 



Formula (26) can be proved as follows. If we define 

u = 8 + X[l - 7(a)], 

then by (21) s = u - X[l - ¥(u)], whence by Burmann's theorem 
(cf. Appendix) for n = 0, 1, • • • we have 









1 



1 - X 



i - *(«; 



(27) 



and the nth derivative can be calculated by using Fan. di Bruno's formula 
(cf. Appendix). On the other hand 



(£?) = l -Xt'(O) = l + xr,, 

\</S/g=0 



and 



(^) = -X 7 ( " +1) (0) = (-l)"Xr n+1 (n = 1,2, ...)■ 

Comparing the above formulas we obtain T„ for every n. 
Finally we remark that, by (25) 



(28) 



(29) 



r 2 = 



Ct-1 



(1 - Xa) 3 ' 



(30) 
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r _ "a _i_ 3a "2~ . v 

3 (1 -Xa) 4 " 1 " (1 -Xa) 6 ' {6U 

T a * - 4- 1QAQ, 2ft3 , l")X 2 a 2 /„ 9 s 

4 (1 -Xa) 5_t " (1 - Xa) 6 ^ (1 - \ a y K6Z) 

IV. THE DISTRIBUTION FUNCTION OF THE DELAY 

Let us denote by r) n the delay of the nth call. If the order of service is 
specified then the distribution function of tj„ is uniquely determined by 
the distribution of (£„ , x»i). If (£n , Xn), n = 1, 2, • • •, is a stationary 
stochastic sequence, then rj„ has the same distribution for every n. In 
the case of the stationary process write P{i?„ ^ x) = W(x) and 

E{e~""'| = Sl(s) 

for each order of service. Define 

W n - [ x n dW(x) (n = 0,1,2, •••)■ (33) 

Iii each case W„ is finite if a H+i is finite. For each order of service 



W -W**5- ™ 



For service in order of arrival 



" J " 3(1 " X«) + 2(1 -\ a y- ' (3o) 

for service in random order 

- 3(1 - X«)(2 - Xa) X (1 - Xa) 2 (2 - Xa) ' k W 

and for service in reverse order of arrival 

jr., - ^ + ^l rs7) 

M - 3d - x«) 2 + 2(i -\ a y {d7) 

(i) Service in Order of Arrival. This case was first investigated by F. 
Pollaczek 4 and A. Y. Khintchine. 5 Cf. also D. V. Lindley/' 

Theorem 2: If \a < 1, if the process is stationary and if service is in 
order of arrival, then the distribution function of the delay of a call is 
given by 

CO 

W(x) = (1 - Xa) £ (\a) k II k *(x) (38) 

fc=0 
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where H k *(x) denotes the k-th iterated convolution of 

- [ [I - H(u)\du ifx ^ 0, 
ifx<0, 



H*(x) = 



(39) 



with itself; H *(x) = 1 if x ^ and H *(x) = if x < 0. 
Proof: Evidently 

80 

W(x) = P + HPjixYHj-Ax) 

where Hj(x) (j = 1, 2, • • •) denotes the jfth iterated convolution of 
H(x) with itself; H (x) = 1 if x ^ and tf„(z) = if x < 0. The 
symbol * denotes convolution. Hence 

fl(*) = P« + t nMMsW-* = Po + m ff )} (4D 

where ?7(s,z) is defined by (7) and P = 1 - Xa. Thus 

1 - \a 



fi(s) = 



! x 1 ~ *(«) ' (42) 



whence (38) follows by inversion. Formula (38) was found by V. E. 

Benes. 

If « n+1 is finite then W« is also finite and is given by 

"•-&r^b r - (43) 

where F„,, is defined by (26). For, 

W n = (-ir(^) s=o (n = 0,l,---) (44) 

and the nth derivative of 0(s) can be obtained by Faa di Bruno's 
formula (cf. Appendix). 

In this case W\ is given by (34), W 2 by (35), and 

_. Xa 4 , xWa , 3X m /k) 

W > = 4(1 - Xa) + (1 - Xa)' + 4(1 - X«)° ' < W 

flemar/c. Let T(x) = W(x)*H(x), i.e., 7'(.t) is the distribution func- 
tion of the sum of the delay and the holding time of a call for a stationary 
process. Define 
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T n = f x n dT(x) (n = 0,1,2, ■■■ ). (46) 

If a n+ i is finite then T„ is also finite and is given by 7\ = W\ + a and 

T n = W n + 5 W„_! (n = 2,3, • • • ) (47) 

A 

where TF„ is defined by (43). Conversely, if we know T, f or j = 1, 2, 
• • • , n then we can obtain W„ by the following formula 



(n = 1, 2, • • • ) . (48) 



_ (-i)-n! r a (-i)Vr ; 

^ n — "r' + S - ji - . 

Formulas (47) and (48) follow from the relationship 

f* e~' x dT(x) = Q(sMs) = (1 ~ x Xa)s + (l - 0(»). (49) 

Finally we also note that the rth binomial moment of the stationary 
distribution of the queue size, i.e., that of {P>} defined by (19), is given 
by 

Br = £(^P 3 = X ^ (r = 0,l, •••). (50) 



For we can easily see that 



i = f e- Xx ^-dT(x), (51) 



P 

'o J! 



whence 



£ r = /" ^ dT(aO = ^ . (52) 

•'o r! r! 

(m) Service in Random Order. The case of exponentially distributed 
holding times was investigated by many authors (cf. Ref. 8), the case 
of constant holding time by P. J. Burke, and the general case by J. F. 
Kingman. 1 The following theorem is due to J. F. Kingman. 10 

Theorem 3: If \a < 1, if the process is stationary and if service is in 
random order, then the distribution function of the delay of a call has 
the following Laplace- Stieltjes transform: 



0(a) = (1 - \a) (l + - f exp|~- f 



dv 



tf(s + X(l - v))_ 



' f ~-%WN ;5:;i 



u — 

u - *(x(: 
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where y(s) is the root with smallest absolute value in z of the equation 

z = *( s + x(l - z)) (54) 

and is given by (22). 

Proof: Under the condition that j(j = 1,2, •■■) calls are waiting in 
the system when a service is about to start, denote by Wj(x) the prob- 
ability that the service of a given call among the j calls starts within 
time x if time is measured from this instant. Define 

0,(8) = r e- sx dWj(x) («(«) ^ 0). (55) 

The distribution functions Wj{x) (j = 1, 2, • ■ •) can be obtained by 
using the following relationships: Wi(x) = 1 and 

Wj(x) = I + (l - -) £ [ f ^" Xu ^ dff(«)l *I^ + *-i(s) (56) 

J \ J/ h=a _ J o Kl 

for j - 2, 3, ■ • • and for .t ^ 0. To prove (56) we take into consideration 
that if the given call will be chosen for service among the j waiting calls, 
then its service starts immediately; if the given call is not chosen for 
service at this time then it must wait during the holding time of the call 
chosen for service, and if during this holding time k new calls arrive, 
then there is an additional delay which has the distribution function 
Wj+k-i(x). Forming the Laplace-Stieltjes transform of (56) we obtain 
the following system of linear equations for the determination of Q,-(s) 
(j = 1,2, ■■•):Sl l (s) = land 

Ms) - 1 + (j - 1) S SW-iU) f c-* +8)z ^f dH(x) (57) 

fc=0 J 1*! 

for j = 2, 3, • • • . The solution of this system is given by J. F. Kingman 10 
in the following form : 

aM ' L ° XP [~ i .-*(.+*(!-.))] ,„, 

» 

du 



u - *(s + X(l - u)) 

where y(s) is defined by (22). By integrating by parts (58) can be 
written in the following equivalent form: 

0,(«) = 1 - (j - 1) I 

7 * d (59) 
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Now the distribution function of the delay is given by 

•In 



w(x) =/nEZ 



/,•! 



dPj(u)\ 
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W j+k (x). (00) 



For, if a call arrives and finds the line free, then its service starts without 
delay. If an arriving call finds j (j = 1, 2, • • • ) calls in the system, then 
its delay is composed of the time needed to complete the current service, 
and if during this time k new calls arrive, then there is an additional 
delay that has the distribution function W j+k (x). Forming the Laplace- 
Stieltjes transform of (('>()) we get 

), (X*)* 



fl( 



*) = Po + i±\ f\- a+s)A ^dPAx)]v J+ M. 

y=l k=ti [_ J Kl J 



(01) 



Putting (58) into (01), we obtain 

a( ' ) - f ' + L CTp [-/.' .-*(.+x(i-.)) ] 

U(s + X(l - u),u) 
'u[u - *(s + X(l - u))] 

where P = 1 — Xa and U(s,z) is defined by (7). Thus 



(02) 



du 



S Jy(s) 

\-f- 



I 



• exp 

u - 1 



dv 



- ¥(s + X(l - v))_ 
u - 1 



(03) 



du 



u - ¥(X(1 - w)) u - ¥(s + X(l - m))_ 

whence (53) follows by integrating by parts and by using the fact that 
y(s) satisfies (54) in z. 

If a„ + i is finite, then W„ is also finite and can be expressed by 



jTW^-c-i)'^). 



for r = 1,2, • • -,n and j = 1, 2, • • • , and by 



\ ds> dz K / s= o, r= i 



(64) 



(05) 



iovj + A: ^ n. By using the following formulas we obtain (34) and (30) : 

fxdWj(x) =^-=^, (00) 

♦'o & — Ka 
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r , AWU v 2a 2 Q- - l)(i - 2) (6 - XoWj - 1) (67) 

i * dF ' (a ° = (2-Xa)(3-2Xa) + (2 - Xa) 2 (3 - 2Xa) » ^ 

and further C7oo = Xa, f/io = Xa 2 /2, C/20 = Xa 3 /3 ; 



\ 2 
t/01 - Xa + 2(1 _ ^ , 

\aj Xa2 , X as 



6 ^2 ' 4(1 - Xa) ' 
and 

X 03 + 3X «2 , X a2 

"02 — -?m ^ — r~ ' 



3(1 - Xa) 2(1 - Xa) 2 

(tit) iSem'ce in Reverse Order of Arrival. The case of exponentially 
distributed service times was investigated by E. Vaulot 11 and the general 
case by J. Riordan 12 and D. M. G. Wishart. 13 Now we shall prove 

Theorem 4: If Xa < I, if the process is stationary and if service is in 
reverse order of arrival, then the distribution function of the delay of a call 
for x ^ 18 given by 

W(x) = (1 - Xa) + X £ e~ u ^L- f [1 - Hi(tt)] du (68) 

where H j(x) denotes the jth iterated convolution of H(x) with itself. 

Proof: Denote by G(x) the probability that the length of a busy 
period in the queueing process considered is <£&. G(x) is given by (23). 
Then we can write that 

W(x) = P + (1 - Po) I} 17V" U$-dH*(u)]*G k (x) (69) 

where G k (x) denotes the fcth iterated convolution of G(x) with itself; 
G (x) = 1 if x ^ and G (x) = if x < 0; H*(x) is defined by (18) 
and P Q = 1 - Xa. For, if an arriving call finds the line free, which has 
probability P , then its service starts without delay; if the line is busy, 
which has probability 1 — P , then its delay is composed of the remain- 
ing holding time of the call being served, which has the distribution 
function H*(x), and if during this time interval k new calls join the 
queue, then there is an additional delay that has the same distribution 
function as the total length of A; independent busy periods. Thus we 
obtain (69). Forming the Laplace-Stieltjes transform of (69) we obtain 
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-(\+a); 



(xxY 



Q(«) = P + ( -l ^2 £ [ 7 ( S )] fc f I-*™ V2JJL [1 //; ., | ,/,. 

a M J o ft! 

— (1 — Xa:j + A — — - 7 - rj 

8 + X[l - 7(8)] 



(70) 



where 7(s) is defined by (22). Since 7(s) satisfies (21) in 2, we get from 
(70) that 



0(g) = (1 - Xa) + 



X[l - y(s)] 



S + X[l - 7 (S)] ' 

By using Lagrange's expansion (cf. Ref. 3, p. 132) we obtain 

X 



(71) 



ft(s) = (1 - Xa) + 



A + s 



(72) 



f (-i)V d H / ggCg + x)]' \ 

ft i! ds'- 1 v (« + x>« y » 

whence (08) follows by inversion. 

If a„+i is finite then IF„ is also finite and we have for n = 2, 3, • • • 
that 

W = T (n ~ 2 + y)! X " F (73) 

"" ft(n-l)l(l-Xa)-^- r "" (73; 

where Y nv is defined by (26). If n =-- 1 then W n is given by (34). For, 
,„ /d n U(s) 



W n = (-!)' 



r/s» 



(n- 0,1,2,.-.). 



(74) 



If we use the notation u — s -\- A[l — y(s)] and s = u — A[l — ¥(«)], 
then we can write that 



fi(s) = (1 - Xa) 4- A 



1 - V(u) 



(75) 



whence by using Burmann's theorem (cf. Appendix) we obtain for 
n = 2, 3, • • • that 



W n = 



(-1)' 



d" 



n — 1 I du* 



1 - X 



1 - ¥(a) 



(76) 



and the nth derivative can be calculated by Faa di Bruno's formula 
(cf. Appendix). 
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In this case W\ is given by (34), W 2 by (37) and 

W - ^ a4 _l_ 3X a-ian . 3X ai (ll\ 

Wi " 4(1 - Xa) 3 + 2(1 - Xa) 4 + 2(1 - Xa) 5 ' k 

Remark. If Po(0 denotes the probability that the line is free at time 
t given that it was free at time t = 0, then we can write that 

W(x) = 1 - [Po(a:) - Po(»)] (78) 

where Po( °° ) = 1 — Xa. 

If G*(x) denotes the probability that the length of a busy period is 
£x for the dual process [H(x),F(x),l], i.e., when the interarrival times 
and holding times are interchanged, then we can write that 

W(x) = 1 - [(?*(») - G*(x)] (79) 

where G*( <*> ) = Xa. 

(iv) Ar Extreme Case. Suppose that in the stationary process the 
service of a particularly chosen call starts when and only when no other 
calls are in the system, i.e., its service is delayed until it becomes the 
only call in the system. Denote by \V*(x) the distribution function of 
the delay of this particular call. 

Theorem 5: 7/ Xa < 1, if the process is stationary, and if a particularly 
chosen call will be served only when no other calls are in the system, then the 
distribution function of the delay of this call is given by TF*(0) = 1 — Xa 
and for x > 

<EHx) = (1 _ Xa)x[1 -G(x)}, (80) 

dx 

where G(x) is the distribution function of the length of a busy period and 
is given by (23). 

Proof: Denote by G n (x) the nth. iterated convolution of G(x) with 
itself. Now we have 



w*(x) = p + E £ 






* G j+k (x). (81 ) 



For, if the particularly chosen call finds the line to be free, then its 
service starts without delay; if it finds j calls (j = 1,2, ■ ■ • ) in the system 
and if during the remaining part of the current service k (k = 0,1, ■ ■ ■) 
more calls arrive, then its delay is composed of the remaining holding 
time of the call being served at its arrival and an additional delay which 
has the same distribution as the sum of j + k mutually independent 
random variables each of which has the same distribution as a busy 
period. 
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Denote by J2*(s) the Laplace-Stieltjes transform of W*(x). By (81) 

00 

a*(«) - P + E n> (*+ x - x 7 (s))[ 7 (s)] y 

r m (82) 

= Po + U(s + X - X 7 (s),y(s)) = Po I 1 - x ! ^ , 

where I\ = 1 - Xa and7(s) is given by (22). We obtain (80) by inver- 
sion. 

If a n+ i is finite then 

W n * = P x n dW*(x) (83) 

is also finite and is given by 

W '* " X0 n +"i r " + ' (' J=1 . 2 .---) (84) 

where r„+i is defined by (25). This follows immediately from (82). In 
particular we have 

w,* = A( * 4 _i_ 5X 2 a 2 g3 , 15XV , . 

4(1 - Xa)<^2(l - Xa) 5 ^4(l - Xa) 6 ' ^ 

APPENDIX 

A.l Burmann's Theorem 

(Cf. Ref. .3, p. 128.) Suppose that the first JV derivatives of f(z) and 
the first N — 1 derivatives of g(z) exist at z = 0. If s = u/g(u) and 
0(0) s* 0, then 

/W=/(0) + E^( d ""'^ W n + <,(,»). (88) 
„=i n! \ o»" ' /„=o 

A. 2 Foa rfi Bruno's Formula 

(Cf. Ref. 14, p. 33.) If 2 = /(v/) where y = g(x), then the nth deriva- 
tive of z = f(g(x)) with respect to x at x = is given by the following 
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Table I 



M 


V 


ji 


it 


h 


3* 


jt 


Cji.j-- in 


1 


1 


1 










Yt. 


2 


1 





1 








Yz 


2 


2 


2 











% 


3 


1 








1 






H 


3 


2 


1 


1 









H 


3 


3 


3 












H 


4 


1 











1 




H 


4 


2 


1 





1 







H 


4 


2 





2 










H 


4 


3 


2 


1 










\i 


4 


4 


4 













Ke 


5 


1 














1 


H 


5 


2 


1 








1 





M 


5 


2 





1 


1 








% 


5 


3 


2 





1 








% 


5 


3 


1 


2 











% 


5 


4 


3 


1 











% 


5 


5 


5 














Mi 



formula 



\ ax" /i=o ^=i \ 



dy" /i/=(/(o) 



where 



Y = 

* n,v 



3l+ZJ2+---+nj n =n 



E 



(89) 



(90) 



provided that the derivatives in question exist. 

Forn g 5 Table I contains all then-tuples (ji ,fa » • ' -»i») satisfying 
the requirements ii + j* 2 + ■ • ■ + J» = » and;'i + 2j 2 + • • • + ftf« = n, 
and in addition the coefficients 



C **--* = jjj,! ■•• j„!(2!)^(3!)« •■• ((n+ l)!)'-» 
which we need in using formula (26). 
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